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Magnetite exhibits unique electronic, magnetic, and structural properties in extreme conditions that are of 
great research interest. Previous studies have suggested a number of transitional models, although the 
nature of magnetite at high pressure remains elusive. We have studied a highly stoichiometric magnetite 
using inelastic X-ray scattering, X-ray diffraction and emission, and Raman spectroscopies in diamond anvil 
cells up to —20 GPa, while complementary electrical conductivity measurements were conducted in a cubic 
anvil cell up to 8.5 GPa. We have observed an elastic softening in the diagonal elastic constants (C u and C 44 ) 
and a hardening in the off-diagonal constant (C 12 ) at ~8 GPa where significant elastic anisotropics in 
longitudinal and transverse acoustic waves occur, especially along the [110] direction. An additional 
vibrational Raman band between the A lg and T 2g modes was also detected at the transition pressure. These 
abnormal elastic and vibrational behaviors of magnetite are attributed to the occurrence of the 
octahedrally- coordinated Fe 2+ -Fe 3+ -Fe 2+ ions charge- ordering along the [110] direction in the inverse 
spinel structure. We propose a new phase diagram of magnetite in which the temperature for the 
metal-insulator and distorted structural transitions decreases with increasing pressure while the 
charge- ordering transition occurs at ~8 GPa and room temperature. 

Magnetite (Fe 3 0 4 ) is the oldest known magnet that exhibits an inverse spinel structure with the chemical 
formula of Fe 3+ ( r ^)(Fe 2+ Fe 3+ )( 0 /i)0 4 where subscripted Td represents the tetrahedral site occupied by 
Fe 3+ ions and subscripted Oh indicates the octahedral site occupied equivalently by Fe 2+ and Fe 3+ ions 1 . 
Due to magnetite's natural geological occurrences in the Earth's crust and mantle as well as technological 
applications in magnetic and electronic materials and nano-composites, numerous experimental and theoretical 
studies have been conducted to understand its electronic, magnetic, and structural properties in extreme pres- 
sures and temperatures (P-T) 2 20 . Of particular interest is the occurrence of the metal-insulator Verwey transition 
at ambient pressure upon cooling below approximately 122 K, the Verwey temperature (TV) 2 , which is strongly 
suppressed with increasing pressure up to approximately 8 GPa 7 " 9,1415 . On the other hand, high-pressure studies 
at room temperature have also suggested a number of electronic, magnetic, and structural transitions in magnetite 
including the inverse-normal spinel transition 91415 , the high-spin to intermediate-spin transition 1619 , the 
enhanced derealization or charge ordering of the 3d electrons of the iron ions 17 ' 20 , the pressure-tuned ideal 
inverse-spinel structure 18 , as well as the occurrence of decreased bond length via changes in the local oxygen 
coordinates 13 . These proposed scenarios differ from each other significantly and entail very fundamental under- 
lying physical mechanisms for understanding the behavior of magnetite at high pressures. 

Based on the analyses of previous high-pressure Mossbauer spectra and refinements of the X-ray diffraction 
patterns 91415 , it has been suggested that magnetite undergoes a transition from the inverse spinel structure to 
the normal spinel structure at pressures of approximately 10-20 GPa and room temperature, in which Fe 2+ 
ions occupy the tetrahedral A site while all Fe 3+ ions occupy the octahedral B site (Fe 2+ (r^)(Fe 3+ Fe 3+ )(o/z)0 4 ). 
The inverse-normal spinel transition would result in a 50% increase in the net magnetic moment from 4 /ig to 6 fi B 
as well as significant local volume changes of the A and B-sites iron ions 91415 . However, this normal spinel 
transition model was not confirmed by many other subsequent studies including Mossbauer and neutron studies 
using highly- stoichiometric magnetite samples 131718 " 20 . On the other hand, combined high-pressure X-ray 
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magnetic circular dichroism (XMCD) and X-ray emission spectro- 
scopy (XES) studies on magnetite have suggested a reduction of the 
net magnetic moment of magnetite by approximately 50% at 12- 
16 GPa, indicating a high-spin to intermediate-spin transition of 
Fe 2+ ions in the octahedral site 16 . However, such a spin-pairing model 
was soon disputed using results obtained from the same XMCD 
technique 19 . 

As clearly pointed out in Verwey's original paper in 1939 and 
many later studies 1 ' 21 " 23 , physical properties of magnetite in extreme 
P-T environments can be very sensitive to the stoichiometry, twin- 
ning, and presence of impurities and defects in the starting sample. 
Specifically, the metal-insulator transition near the TV would occur 
continuously or discontinuously depending on the stoichiometry of 
the magnetite samples 1>21 ' 22 ; for a chemical formula Fe 3 ( 1 _5)0 4 where 
5 is the degree of non- stoichiometry, the Verwey transition is of 
first order or second order when the absolute 5 value is less than 
or more than 0.004, respectively. Twinning domains, on the other 
hand, have hampered the diffraction studies of the crystal struc- 
ture 23,24 . Thus, care needs to be taken in preparing and characterizing 
the stoichiometry and avoiding twinning of the starting sample 
in order to better understand the behavior of magnetite at extreme 
environments. It should be noted that natural magnetite samples 
are intrinsically highly non-stoichiometric and can develop 
twinning domains 1 . Unfortunately, natural or synthetic magnetite 
samples have been commonly used in many of the literature studies 
without characterizing and reporting their non-stoichiometric 5 
values as well as the degree of twinning domains 1 , making a coherent 
comparison and understanding of the literature results extremely 
difficult. 

A survey of previous studies at high pressures using highly stoi- 
chiometric starting samples shows that magnetite undergoes a trans- 
ition at approximately 7-8 GPa and room temperature 1318 ' 20 . The 
abnormal changes in the hyperfine interaction parameters at 
7 GPa have been interpreted as the occurrence of the decreased local 
oxygen internal coordinates (the decreased Fe-O bond length) at 
high pressures and room temperature 7 ' 813 . However, high-pressure 
neutron diffraction results and ab initio calculations do not reveal 
any detectable change in the oxygen atomic parameter, yet show 
significant but continuous decrease of the magnetic moments at both 
the A and B sites to at least 10 GPa 17 . These results, nevertheless, rule 
out the possibility of a pressure-induced transition to the normal 
spinel structure at high pressures, and point to the possible occur- 
rence of the charge-ordering of the iron ions at the B site 16 . 
Furthermore, a change in the curvature of the electrical conductivity 
at 6 GPa and room temperature has been attributed to be a result 
of the occurrence of the pressure-tuned ideal inverse-spinel 
structure 18 . 

To decipher the behavior of magnetite at high pressures, here we 
have studied a highly stoichiometric and twinning-free single-crystal 
magnetite sample with a chemical formula Fe 3 ( 1 _5)0 4 , where |8| = 
0.00018 (Fig. SI), using high-energy resolution inelastic X-ray scat- 
tering (HERIX) 26 ' 27 , X-ray diffraction (XRD), Raman spectroscopy, 
and X-ray emission spectroscopy (XES) 28,29 measurements in dia- 
mond anvil cells (DACs) up to approximately 20 GPa (See 
Methods and SI for details) 22 ' 25,30 . Complimentary electrical conduc- 
tivity (EC) measurements were also conducted in a cubic anvil cell up 
to 8.5 GPa at room temperature (See SI for details). The HERIX and 
XRD results were used to derive full elastic constants of the single- 
crystal magnetite, while the Raman, XES, and EC results were used to 
understand the vibrational and electronic properties of the sample. 
These results show abnormal acoustic phonon and optical vibra- 
tional behaviors occurring at approximately 8 GPa. Using these 
experimental results and literature data, here we attribute the origin 
of the abnormal elastic and vibrational properties of magnetite at 
approximately 8 GPa to the local site distortions associated with 
charge-ordering in the inverse spinel lattice. 



Results 

High-pressure HERIX spectra of the single -crystal magnetite platelet 
with the pre-oriented (100) crystallographic plane were measured up 
to 19 GPa in a DAC for the lower momentum transfer (Q) part of the 
longitudinal acoustic (LA) and transverse acoustic (TA) phonon 
branches along [100] and [110] directions (See Methods and SI for 
details; Figs. S2-4) 26 ' 32 . X-ray diffraction patterns were also collected 
to determine the crystallographic orientation matrix and to calculate 
the lattice parameter and the density of the sample as a function of 
pressure (Figs. S5; Table SI). The XRD patterns showed sharp reflec- 
tions of the cubic inverse spinel structure without any observable 
structural transition, twinning domains, or volume collapse up to 
20 GPa within the experimental uncertainties (Fig. S5). Using a lin- 
ear fitting of the measured inelastic energy shift (E) as a function of Q 
near (800) and (440) reflection spots 26 , we have determined full 
elastic constants (C n , C 12 , and C 44 ) along principle [100] and 
[110] axes for the longitudinal (Vp) and transverse (Vs) acoustic 
waves based on the following equations: 

^[100] = (C77/ P ) 1/2 

Vs m = (C 44 /p) 1/2 

Vp [U0] = [(C n + C 12 + 2C 44 )/2p} 1/2 

Vs [110] = [(C 77 -C 72 )/2p] 1/2 

where p is the density and the subscript [uvw] indicates the crystal- 
lographic direction of the acoustic phonon wave propagation 
(Fig. I) 32 . Using these elastic constants, we have also derived adia- 
batic bulk modulus (K s ), shear modulus (G), and aggregate Vp and 
Vs of the single crystal 30 . 

The C n and C 12 elastic constants show monotonic increases with 
increasing pressure up to approximately 7 GPa, while the C 44 
decreases with increasing pressure (Fig. 1). Between 7 and 8.5 GPa, 
abnormal changes in the constants and sound velocities occur in 
which the diagonal elastic constants C n and C 44 drop by approxi- 
mately 6-7%, respectively, while the off-diagonal elastic constant C 12 
jumps by 17% (Fig. 1). The aggregate Vp, Vs and G also decrease 
significantly within this pressure interval. However, the K s continu- 
ously increases with increasing pressure normally up to 19 GPa and 
does not show any visible change at the transition pressure range for 
other elastic constants; this behavior of the K s as a function of pres- 
sure is consistent with the compressional pressure-volume curve 
derived from the X-ray diffraction spectra as a function of pressure 
(Fig. S5). Between 16 GPa and 20 GPa, the C 12f G, and Vs also 
significantly decrease with increasing pressure. The softening of 
these constants between 16 GPa and 20 GPa in the inverse spinel 
phase can be attributed to its approaching the previously reported 
orthorhombic phase at approximately 22 GPa 3 ' 5 ' 6 . 

To understand the abnormal behavior of the elasticity as a func- 
tion of pressure in detail, we have also calculated the Vp, Vs l5 and Vs 2 
as a function of the crystallographic direction at high pressures 
(Figs. 2, S4) 32 . The Vp and Vsi exhibit significant changes especially 
along the [110] direction between 7 and 8.5 GPa with a significantly 
increased anisotropy factor (A) (Fig. 2). The V s splitting anisotropy 
factor (A Vs ) jumps from a few percent to 17% while the V p anisotropy 
factor (Ayp) increases from approximately 1% to 4% (Fig. S4), indi- 
cating that the single-crystal magnetite becomes highly elastically 
anisotropic at pressures above approximately 8 GPa. We note that 
X-ray diffraction patterns of the sample in the soft Ne medium 
showed sharp and round diffraction spots with well converged lattice 
parameters, ruling out the possibility that the observed anisotropy 
was induced by the non-hydrostaticity in the sample chamber. 
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Figure 1 | Elastic constants and sound velocities of single-crystal magnetite at high pressures. Solid circles: this study; open circles: Reichmann and 
Jacobsen (2004) 12 . Black lines represent simple polynomial fits to the experimental data, and are only meant to guide the eyes. Vertical lines indicate the 
transition pressure range at approximately 8 GPa. Vp and Vs represent the aggregate compressional and shear wave velocity, respectively, calculated from 
the single- crystal elasticity data. 



To decipher the observed abnormal transition in the elasticity of 
magnetite, we have also conducted complimentary Raman spectro- 
scopy, synchrotron Fe Kp XES, and the bulk electrical conductivity 
(EC) experiments at high pressures (See Methods and SI, Fig. 3, S6- 
7) 28 ' 29 . Analyses of the high-pressure Raman spectra show three dis- 
tinct Raman bands that can be assigned to A lg) T 2g , and E g modes of 
the cubic inverse spinel structure up to 21 GPa (Fig. 3) 10 > 11 > 33 > 34 . 
Raman frequencies of these modes increase almost linearly with 
increasing pressure up to the highest pressure of 21 GPa, although 
the relative intensity of the T 2g mode with respect to the A lg mode 
starts to increase with increasing pressure at approximately 8 GPa 
(Fig. 3a). Most significantly, an additional Raman band occurs in the 
spectral range between the A lg and T 2g bands, and its Raman fre- 
quency also increases linearly with increasing pressure, similar to 
that of other magnetite Raman bands (Fig. 3). 

Analyses of the synchrotron Fe Kp X-ray emission spectra using 
the absolute and relative integrated areas methods proposed 
previously 28,29 show that the total spin momentum (S) of all iron 
ions slightly decreases with increasing pressure without any 
observable discontinuous changes up to 19 GPa (Fig. S6). Based 
on these spectral analyses, the slight continuous decrease in the 
intensity of the satellite Kp can be attributed to the spectral broad- 
ening effect, instead of any spin-pairing transition of iron in mag- 
netite (Fig. S6) 1619 . Electrical resistance of the single-crystal 
magnetite as a function of pressure measured in a cubic anvil cell 
continuously decreases with increasing pressure (Fig. S7), and a con- 
tinuous change in the curvature can be identified at approximately 
7 GPa, consistent with previous results 7 ' 818 . The electrical resistance 
results indicate that magnetite becomes a better metallic conductor at 
higher pressures. 



Discussion 

Previous high-pressure ultrasonic interferometry and Raman mea- 
surements 10 " 12 on highly non-stoichiometric natural magnetite sam- 
ples did not reveal any changes in elasticity nor the occurrence of the 
additional Raman band at the pressure range of 7-8 GPa where we 
observed the occurrence of the abnormal elastic and vibrational 
behaviors on a highly stoichiometric magnetite sample. Our elasticity 
results also differ from the high-pressure ultrasonic interferometry 
results 12 on a natural single- crystal magnetite by approximately 10%, 
although the differences become smaller at higher pressures (Fig. 1). 
Thus, it is evident that the differences in these observations can be 
attributed to the stoichiometry of the samples used in these studies; 
while our sample was highly stoichiometric with a chemical com- 
position of Fe 3 ( 1 _5)0 4 where |5| = 0.00018, 25 the natural magnetite 
sample used in the ultrasonic study was likely non-stoichiometric 12 . 
Together with the inconsistency in the literature results on other 
physical properties of magnetite at high pressures, it is clear that 
the high-pressure behavior of magnetite is highly sensitive to the 
stoichiometry of the starting sample and that the abnormal transition 
at approximately 8 GPa only occurs in highly- stoichiometric sam- 
ples. Therefore, experimental results obtained from highly stoichi- 
ometric samples are of direct relevance to understanding the high- 
pressure behavior of the sample. 

Based on our X-ray emission results, we can clearly rule out the 
spin -pairing transition model of the iron ions in magnetite at high 
pressures 1619 . The inverse to normal spinel transition model 91415 is 
also unlikely to contribute to our observed abnormal elastic and 
vibrational behavior, because such a transition would result in sig- 
nificant changes in the local site symmetry and volumes of the lattice 
as well as in the electrical transport property; all of which were also 
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Figure 2 | Compressional and shear wave anisotropics of single-crystal magnetite at high pressures. The notation for the two polarized shear waves is 
defined as V Si being slower in velocity than V S2 . Principle axis [100] is along the horizontal line while the vertical direction is [010]. Vp and V s 
velocities exhibit an enhanced anisotropy starting at 8.5 GPa. Specifically, the V S1 and Vp change dramatically along the [110] direction across the 
transition at approximately 8 GPa. 



not observed here nor in previous studies using highly stoichiometric 
samples 7 ' 8 ' 13 ' 17 ' 18 ' 20 . 

Based on the group theory and lattice-vibration analyses, the 
inverse spinel structure of magnetite in the primitive unit cell con- 
tains Fe and O atoms in three distinct site symmetries (43m (T d ), 3m 
(D 3d ) and 3m (C 3V )), in which the T d and C 3V sites occupied by Fe 3+ 




Raman Shift (cm" 1 ) Pressure (GPa) 



Figure 3 | Vibrational Raman modes of single-crystal Fe 3 0 4 at high 
pressures, (a) Representative Raman spectra; (b) derived Raman 
frequencies. A new Raman band occurs between A lg and T 2g modes 
(orange line) starting at approximately 8 GPa. 



and 0 2+ ions, respectively, contribute to the five Raman- active 
modes of the Fe0 4 tetrahedron including one A lg , one E g , and three 
T 2g (Fig. 3) 1 0' 11 33 ) 34 The A ^ mode at _ 668 cm -i arises from the 

symmetric stretch of the O atoms along Fe-O bonds, the T 2g mode 
at —538 cm" 1 is the asymmetric stretching of the Fe and O atoms, 
and the Eg mode at —306 cm" 1 comes from the symmetric bonds 
between Fe and O atoms. Due to the delocalized nature of the indis- 
tinguishable Fe 2+ and Fe 3+ 3d electrons in magnetite at ambient 
conditions, the octahedral site retains the local 3m (D 3d ) symmetry, 
which would otherwise exhibit a lower symmetry due to the occu- 
pation of both Fe 2+ and Fe 3+ in the more "ideal" inverse spinel 
structure 33 ' 34 . Such charge transfer also contributes to the "half 
metal" electrical conduction behavior of magnetite through the 
overlapping integrals between Fe 2+ and Fe 3+ ions via the octahedral 
network above the Verwey transition temperature at ambient 
pressure 7,818 . 

Recently, it has been reported from detailed X-ray diffraction 
structural analyses of a single-crystal magnetite that Fe 3+ -Fe 2+ - 
Fe 3+ ions at the B-site form linear three-Fe-site units, called trimer- 
ons, with the charge order and three-site distortions along the [110] 
direction below TV (Ref. 23). The formation of the trimerons can 
induce substantial atomic displacements and affect the coupling 
between the large electrical polarization and the magnetization 23 . It 
has also been proposed that these trimerons may become impor- 
tant quasiparticles in magnetite above the Verwey temperature 23 . 
Examinations of the measured Vp and V s velocities as well as the 
derived Cij of the single- crystal magnetite indeed show much more 
dramatic changes in the V S i and Vp velocities along the [110] dir- 
ection across the transition at approximately 8 GPa than other direc- 
tions (Figs. 1-2). The pressure-induced isostructural transition 
at 8 GPa with enhanced anisotropic elasticity suggests that the 
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Figure 4 | Phase diagram of Fe 3 0 4 at extreme pressure and temperature. 

The diagram is constructed using our results as well as literature 
results 2 ' 4,7 " 913 " 1518 ' 20 ' 21 . Magnetite undergoes the Verwey transition (Tv) at 
ambient pressure where structural, charge ordering, and metal-insulator 
transitions occur. The temperature for the metal-insulator transition 
(solid grey line) 7 ' 8 and the distorted structural transition (solid black line) 14 
decreases with increasing pressure, while the temperature for the charge 
ordering transition (blue dashed line) increases with increasing 
pressure 9 ' 13,18,20 . 

Fe 3+ -Fe 2+ charges at the B-site become ordered at pressures above 
8 GPa. It should be noted here, however, that the Verway transition 
is accompanied with the trimeron ordering, the monoclinic struc- 
tural transition, and the metal-insulator transition at ambient pres- 
sure 1 (Fig. 4). As the transition to the trimeron ordered phase is 
separated from the metal-insulator transition at increasing pressures 
(Fig. 4), there is no longer a dramatic change in the electrical trans- 
port property. Ordered trimerons would result in a less symmetrical 
yet "ideal" inverse spinel structure with a small displacement (dis- 
tortion) along the [110] direction which would contribute to addi- 
tional Raman modes and would affect the elastic constants and 
sound velocities along the direction more significantly than other 
directions 33 . Thus, the iso- symmetric electronic transition model 
from the charge delocalized state to the ordered trimeron in the 
octahedral site of magnetite at approximately 8 GPa is most consist- 
ent with our observations and literature results from highly stoichi- 
ometric samples. Our observed transitional pressure is also 
consistent with the Mossbauer studies on highly stoichiometric mag- 
netite samples, which showed significant changes in the hyperfine 
parameters, especially in the octahedral sites 913 " 15 . 

The strength of the electron -phonon interaction in magnetite 
plays an important role in its unique electronic and magnetic prop- 
erties at ambient conditions. It is conceivable that the occurrence of 
the ordered electronic transition at approximately 8 GPa signifi- 
cantly affects the electron-phonon interaction, resulting in the 
abnormal acoustic phonon behavior in the inverse spinel structure 
observed in the HERIX measurements. The softened C u (diagonal 
compressional component) and C 44 (diagonal shear component) and 
enhanced C 12 (off-diagonal component) observed here would thus 
hold the key to deciphering this localized electronic transition. Since 
previous studies have suggested that the transition in magnetite at 
approximately 8 GPa and 300 K can extend very close to the Verwey 
transition at approximately 122 K and ambient pressure 7 " 91415 
(Fig. 4), future investigations at high pressure and low- temperature 
conditions are needed to decipher the relationship between the 
Verwey transition and the transition to the trimeron ordered phase 
in the inverse spinel structure. 



Methods 

Synthesized single-crystal magnetite with a highly stoichiometric chemical com- 
position of Fe 3 ( 1 _§)0 4 , where |5| = 0.00018, was used for all of our high-pressure 
experiments (SI; Figs. Sl-7) 22 ' 25 . A millimeter- size piece of the single-crystal mag- 
netite was oriented in the (100) plane using single-crystal X-ray diffraction (XRD) at 
the Texas Materials Institute of the University of Texas at Austin. The oriented 
sample of approximately 1 50 um in diameter was then double- side polished to 40 um 
in thickness. A Re gasket of 5 mm in diameter and 250 um in thickness was pre- 
indented to 1 5 GPa using a pair of diamond anvils having 500 um flat culets in a short 
symmetric DAC. A hole of 250 um was drilled on the pre-indented area of 
approximately 75 um thick, and then used as the sample chamber. The sample was 
loaded into the sample chamber, together with two ruby spheres as the pressure 
calibrant 29 . Ne gas was loaded into the sample chamber as the pressure medium for 
inelastic and X-ray diffraction experiments using the Gas Loading System at the 
Mineral Physics Laboratory of the University of Texas at Austin. 

High-energy resolution inelastic X-ray scattering (HERIX) 26 ' 27 and XRD experi- 
ments were conducted in a short symmetric diamond anvil cell (DAC) up to 19 GPa 
at 300 K at Sector 30 of the Advanced Photon Source (APS), Argonne National 
Laboratory. The DAC has an opening of approximately ±40° for the experiments in 
the axial direction. An incident X-ray beam with an energy of approximately 
23.724 keV and an energy bandwidth of 1.5 meV was focused to a beamsize of 15 X 
35 um at the sample 27 . An online image plate diffraction detector was used to orient 
the crystal with the (800) or (440) diffraction spot aligned in the horizontal plane. The 
scattered inelastic signals of the sample were collected and analyzed using a spher- 
ically-bent silicon crystal analyzer of the (12 12 12) reflection and a detector working 
very close to back reflection (89.98°) in the energy range of ± 10-12 meV with a step 
size of 0.5 meV 31 . Typically, the HERIX spectra for the longitudinal and transverse 
acoustic phonons were collected along the [100] and [110] directions with a spectral 
energy resolution of approximately 1.5 meV and a momentum resolution of 
0.2 nm" 1 (Figs. S2, S3). The collection time for each energy scan was approximately 
2 hours and about 3-4 spectra were co-added for each set of experimental conditions. 
XRD diffraction spots of the sample were also collected and used to determine the 
sample densities as a function of pressure, while an online ruby system was used to 
measure the ruby fluorescence for pressure determinations (Fig. S5). 

High-pressure Raman spectroscopic measurements on a (100) -oriented single- 
crystal magnetite were conducted in a symmetric DAC having a pair of 200 um 
diamond anvils which consisted of a natural ultralow- fluorescence diamond anvil 
and a synthetic ultralow- fluorescence diamond anvil produced by Sumitomo 
Corporation. Since the Raman signals of the opaque magnetite were extremely weak, 
the use of the synthetic diamond with extremely low fluorescence background sig- 
nificantly enhanced the signal-to-noise ratio of the Raman spectra. Similar to the 
HERIX experiments, the sample was loaded into the sample chamber of a DAC with 
two ruby spheres as the pressure calibrant 30 and Ne gas as the pressure medium. A 
Renishaw in Via Raman Microscope system equipped with a 532 nm laser and a long- 
working distance lens of 20 X was used for collecting the Raman spectra with a typical 
data collection time of 30 mins. The HERIX and Raman spectra were analyzed using 
the OriginPro 8 Program. 
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